arXiv:1505.00868vl [math-ph] 5 May 2015 


Formation of compression waves 
with multiscale asymptotics 
in the Burgers and KdV models 

S.V. Zakharov 

Institute of Mathematics and Mechanics, 

Ural Branch of the Russian Academy of Sciences, 

16, S.Kovalevskaja street, 620990, Ekaterinburg, Russia 


Abstract. The Cauchy problem for the Burgers equation with a small dissipation and 
an initial weak discontinuity and the Cauchy problem with a large initial gradient for a 
quasilinear parabolic equation and for the Korteweg-de Vries (KdV) equation are considered. 
Multiscale asymptotics of solutions corresponding to shock waves are constructed. Some 
results can also be applied to rarefaction waves. 

Keywords: Burgers equation, KdV equation, Cauchy problem, multiscale asymptotics, 
weak discontinuity, shock waves, renormalization. 

Mathematics Subject Classification: 35K15, 35K59, 35Q53. 


1 Wave with an initial weak discontinuity 


A simplest model of the motion of continuum, which takes into account nonlinear effects and 
dissipation, is the equation of nonlinear diffusion 


du du d 2 u 
dt U dx £ dx 2 ’ 


£ > 0 , 


( 1 . 1 ) 


for the first time presented by J. Burgers [lj. This equation is used in studying the evolution 
of a wide class of physical systems and probabilistic processes, for example, acoustic waves 
in fluid and gas M- The problem of shock wave formation (for equation (II.Tj) with e > 0 ) 
is investigated in [3, UEj- A perturbed weak discontinuity is considered in [2]. The problem 
of shock wave propagation for the Hopf type equation (e = 0) is studied in [71 [S]. 

Below, we briefly state results of paper [4] for the following initial data: 


u(x, — 1 , e) = — (x + ax 2 )0(—x), 16 R 1 , ( 1 . 2 ) 

where a > 0, Q(x)\ x ^o = 1, @(x)| x <o = 0. The solution of the limit (e = 0) equation is found 
using the method of characteristics: 


u 0 (x,t) = 


2a(l + t)x + t + \Jt 2 — 4a(l + t)x 


0, 


2 a(l + 1) 2 


, x < s(t) 0 (t), 
x > s(t)Q(t). 


The function uq has a jump discontinuity on the line x — s(t), where 


. , (3 1 + 4) 3//2 - 9t - 8 

s(i) =-MTT7)-' 


s(t ) = 


3 1 2 


271 3 


16 a 128a 


+ 0(f 4 ), t-> 0. 


(1.3) 
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The solution of problem fll.lH - fll.2H is given by the following expression: 


u(x,t,e) = -2e[^(x,t,e)] 1 '4/ X (x,t,£), 


(1.4) 


where 


V(x,t,e) = 


2 yj 7r(l + t) 


exp < - 


(■ x~y) 2 v 2 ay 


+ V + 


e l 4 (1 + t) 4 6 

(x - y ) 2 


2a/7t(1 + t) J ^ L 4(l+f) 
v o 

Taking into account these formulas, we write expression (11.41) in the form 

T°(a;, t, £) 


1 


dy + 


dy. 


u(x, t, £) = 


T (x, t, £) + T + (x, t, £) 


where 


T (x, t,£) — / exp 


F ( y,x,t ) 


(1.5) 


dy, F (: y,x,t) = 


(: y-xr , IJ , ay 


+ V + 


4(1 + t) 4 6 ’ 


T + (a;, t, £) = 2y/£7r(l + t) — + — exp 


x 


e 


2 y/T+i 


4(1 + t) 


+ y/efTTi) / — z~dz, (1.6) 


C = x/y/e, 


4>°(x,t,£) = / (y Fay 2 ) ex p 


F (y,a;,t) 




The integral T (x, t, e) can be written in the form of the sum 

y + (x,t) o 

ty~(x,t,£) = Vj(x,t,e) + %(x,t,£) = J e F/£ dy+ J e F/s dy , 

—°o y+(x,t) 


where y + (x,t) is the point of a local minimum of F ( y,x,t ) in the variable y. The asymp¬ 
totics of 4>j(x,t, e) is found by Laplace’s method. Therefore, we need expressions 


V T (x,t ) 


t ± y/t 2 — 4ax(l + t) 
2a(l + 1) 


(1.7) 


for extremal points of the function F~ in the variable y (for fixed x and t) and values of the 
function F~ and its derivatives at these points. 

Substituting y = y T (x,t ) into expressions for F~(y,x,t) and F~(y,x,t), we have 

F yy>x,t) = T’ where = “ 4ax(l + t), (1.8) 
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F {y T (x,t),x,t) = 


Thus, 


where 


(x, t, e) = \Z2 he 


H(x, t ) = 


exp 


xt + 3ax 2 (l + t) - [■ R ( a; » t )] 2 

12a(l + 1) 2 ' 

F~(y~(x,t),x,t) 


(1.9) 


\J~ F yy(y {x,t),X,t) 

T ■ 1 , . X — S(t) 

Let us introduce a stretched variable a = - 


'Y^e?'y j a 23 [H{x,t)]*i +l , (1.10) 

i=o 

7o = 1, 7i = 5/24. 


, where s(t) is from (11.3|) . We make 


the change x = s(t ) + ea in expressions for R(x,t), H(x,t ), and y ± (x,t ): 

r / ,\i2 N — l 

R(x,t ) = [ -^ + T(ae) m r m (t) + 0(e aN ), r m (t) = O^ 1 " 2 -), (1.11) 

Z= S (t)+£<7 -< ^' 


7TI— 1 


F(:r, £) 


2(1 + 1 ) 


N-l 


x=s(t)+ecr V l?(s(t) T £<7, t) 

where r(£) = ^/3t + 4 — 2\/3 1 + Z, Co = 1, ci = 9; 

IV -1 


+0(0. (1.12) 


m=0 


[r(t)} 


y ± (x,t ) 


x=s(t)+£cr 


+ o(£“ w ). i™(«) = 07 


—2 m\ 


(1.13) 


m=0 


Using expansion (11.111) . let us pass to variable a in formula (11.91) 


N—l 


F {y (x,t),x,t) 


F (y + (x,t),x,t) 


X=s(t)+Ea 


—y(t)<JE + ^(£a) m F_ m (t) + 0(t 3 £ aN ), 


1 [r(t)] 


m =2 

6 N-l 


x=s(t)+ecr 324a 2 (1 T t)^ 

v ' m= 1 

where F^Jt) = 0(t 3 ~ 2m ), p(t) = s'(() = 


V*) 


(3t + l)\/3tT4 - 2 


(1.14) 

(1.15) 


18a(l + 1)' 2 


Using expansions (11.101) and (I1.12j) . we obtain 


T;(x,t,e) = V2 7T e exp 


F (j/ (x,t),x,t) 


a/ 6(1 +t) 

r(t) 


iV-l 


1 + t) + 0(e 


ck/Va 


ra=l 


where G II m . To describe the behavior of coefficients of some expansion, we used the 
class of functions 

n m = {T G C°° : \*(a,t)\^P m (<Tt- 2 ,t- 3 ) (cr,t) Gl'x (0, T)} 

(. P m is a homogeneous polynomial of degree m). 

Now, let us consider ^f^(x,t,e). For t > e 1//4_A 

|T"(x,t,£)| ^ \y + (x,t)\ exp [F“(0,x,t)/e] = O (exp[-e _A °]) , A 0 > 0. 
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For t < £+ 4 a we represent it in the form sum: 

- g v 


^ b (x,t,£)= / e F ,£ dy + / e F /£ dy, 


y+(x,t) —9^ 

where g > 0, 3/2 < 7 < 3/(1 — a) — 2. The first integral is estimated as follows: 

-gt-y 


e F t e dy 




< \y + (x,t )| exp [F (-gF,x,t)/£\ = O (e 9lf2+7/e ) 


The inequality 7 < 3/(1 — a) — 2 provides an exponential smallness of the remainder in the 
domain 

= {|cr| < f 2 £ 4 “ _1 , £ 1_a < t 3 < const, 0 < a < 1 }. 

From the inequality 7 > 3/2 we conclude that for — gt 1 ^ y ^ 0 and t < £+ 4_A 


5(y,t) = - 


ty 2 ay 


+ 


£ [4(1 + t) 6 J 


= O (e Al ) , where Ai > 0 . 


This allows one to expand e F + in the second integral into a Taylor series in the small 
parameter S(y,t ): 


T b = exp 


c 2 


4(1 + *) 


N~ 1 id 71 - 1 )/ 2 ] 1 + \m+l 

E-’” JLt? +<?++ 


X 


|^m=l 1=0 

where, in particular, a 1 ,0 = 2, 02,0 = 4. Let us pass to variable a: 


4b = exp 


c 2 


2e{l + t) 


N -1 


X 


4(1 + *) 

Thus, we obtain 

= v*27T£exp 


) + 0(e 

m =2 


aN \ 


e n m . 


F (y (x,t),a;,t) 


a/6(1 +1) 

r(t) 


2 s (1 + f) 




exp 


c 2 


4(1 + *) 


exp 


c 2 


4(1 + *) 


O ( t £ 2 °) 


(1.16) 


The function 4+ can be represented in the form of the sum: 4+ = 4+ + 4/(, where 


y + (x,t) 


4 '°(x,t,e) = / (y + ay 2 )e F /e dy, ^° b (x,t,£) = / (y + ay 2 )e F ,£ dy. 


y+(x,t) 


Functions 4+ and T b are studied in exactly the same way as 4/ s and 4/ b , respectively. Thus, 


4/ b (x, t, e) = exp 


c 


4(1 + *) 


'JV-l 


E^U^J + ob 


(y.N\ 


m =2 


e n m, 
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\l/°(x,t, e) = —\Z2ite 


exp 


F (y (x,t),x,t) 


j =o 


where 


° t + R(x,t) + 2ax(l + t) 

j 2a(l + 1) 2 { h 


*?(»■ t) =«[*(», t)F-> +«)i w -+ 7 +1+f) [*(», or 1 - 


Passing to variable a, we obtain 


where G n m . 


a/6(1 +£) 

r(t) 


N -1 


2/ J (<) + E £ ’”<»(‘ r '‘) + 0 ^ 


aiV\ 


m=l 


(1.17)' 


Using expression ('ll.61) . (11.161) . and (11.171) . we arrive at the following result [5j. 
Theorem 1. In the domain 

Hq = {|cr| < £ 2 £ 4 “ -1 , E l ~ a < t 3 < const, 0 < a < 1} 

for the solution of problem (11.11) (11.21) , there holds the asymptotic formula 

N-l 


u 


( x, t,e) = ^2 £ p/2 h p (a, C, t) + 0(e N ), 

p =o 


(1.18) 


where N —y oo as N —y oo. The leading terms of this expansion is 

2/4*) 


h 0 (a,C,t) = 


1 + K(t) exp(/i(t)cr) 


-C/(2 v / l+i) 


1 + 


2 ^ 


-dz 


where K(t) = ^1 + ^ + 0(t 2 ) ) . 


Expansion (j 1.18 j) shows the presence of an additional scale in the form of the stretched 
variable ( = x/\fe. 


2 Large initial gradient 

In this section, we consider the Cauchy problem for a more general quasilinear parabolic 
equation mmm-- 


du dipiu) d 2 u 

r, + o — 2 ’ £>0, 

at ox ox 2 

(2.1) 

u(x, 0, e, p) = z/(xp _1 ), igR, p > 0. 

(2.2) 
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We assume that the function ip is infinitely differentiable and its second derivative is strictly 
positive. The initial function v is bounded and smooth. The interest to the problem under 
consideration is explained by applications to studying processes of shock waves formation [31 

HUGE!. 

The asymptotics of the solution as e —» 0 and /j = p/e —> 0 in the leading approximation 
has the form 


u(x,t,e,p) 




+ r (?t) + °b I/2 ‘‘id- 


where ho, Ro,o,o and T are known functions. Here, a multiscale character of asymptotics 
arises from the beginning because of another small parameter p. 

The behavior of the solution of problem fl2.1| ~ fl2.2D is mainly determined by the solution 
of the limit problem 


du d<p(u) . . ii/ n 

— + = 0 , u(x, 0 ) = J 0 


v 0 , x < 0, 
Uq, x ^ 0. 


dt dx 

For Vq > Vq using the method of characteristics, we find its generalized solution 

<^Oo + ) - 


Uo t o(x, t) — 


u 0 , X < ct, 
Uq, X > ct, 


c = 


v o -Vq 


This solution is discontinuous on the line of the shock wave x = ct. 
In paper HU, for problem (j2.ip ~ fl2.2p outer expansions 


oo m— 1 


U + (x,t,£,p) = V+ + 


m= 1 n =0 
oo m —1 


U (X, t,E,p)=U o + ^ pm n£nU rn,n( X i *) 

m=l n=0 

are constructed in domains 

Qq — {(x,t) : x > ct + £ 1 ~ s °, 0 < 5 0 < 1} 

and 


Q 0 = {(x, t) : x < ct — e 1 <5 °}, 


respectively, where 


a^ ns are constants. 


m —1 




x = 


t s 

m.n.s 




(2,3) 


(2.4) 

(2.5) 


Now let us construct an asymptotic solution in a neighborhood of the line of discontinuity. 
First, rewriting outer expansions (12.4|) and (12. 5 ji in terms of the inner variable 


x — ct 


a = 


£ 


we obtain 


In/ 2] 


C/*(d + ea,t, e, p) = i4 + Yi Yi np t- 2 m( a ). 


( 2 . 6 ) 


n= 1 m= 1 
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where Pn- 2 m ( a ) are polynomials of degree n — 2m. Taking into account the structure of 
series ( 12 . 6 p and the structure of the inner expansion 

OO 

H = ^2/d n h n (x/p,uj), a; = te/p 2 , 

n =0 


oo m 

.-m/2 


Icr| + U) —> OO, 


K = u n ' 2 Y u- m/2 u) lR n, m ,l ( Y~r 

m =0 1=0 ' 

we will construct an asymptotic solution of equation ( 12 .ip in the domain 

Q 3 = {(x,t) : \x - ct\ < e l ~ 6z , t > £ 1-73 , 0 < y 3 , 5 0 < S 3 < 1} 


in the form of the series 

oo n 

V(a, t, p, e) = v 0 (a) + Y Y /x m £ n ~ m Y (In fi) p (In e) q v m}n - m! p !q (a, t), 

n= 1 m =1 O^p-hq^n 

where 

v 0 (cr) — A (cr + x) (2.7) 

is the leading term, the function A is defined by the formula 


/ 


K + +q> )/ 2 


dv 

<p(v) — cv — b 


= <r, 


c = 


^o + - 


6 = 


) - ^0 7>0( 




^0 - ^0 


x is constant, which should be determined by the matching procedure. 
In a standard way we arrive at the system of equations 


where 


<9 2 u 0 5<p(w 0 ) 

<9a 2 C 5a 5a 

(2.8) 

Lz v l,0, p,q — 0, 

(2.9) 

1 _ m,p,q . dQm,n,p,q 

l -‘3 v m,n-m,p,q ~ r^_ d(J ’ 

(2.10) 

5 2 u 5u dip'(vo)v 

3W “ 5a 2 + C da da ’ 

(2.11) 

y ^ (r) (u 0 ) V- TT 

Wm,n,p,q | / , || Ur mk,n k -m k ,p k ,q k ; 

(2.12) 


r=2 6 k =1 


© = {( m k ,n k ,p k ,q k ) : Y m k = m ,Y nk = n,YP k = P, ^ ( lk = q}- 

k= 1 /c=l /c=l /c=l 

We should hnd solutions of the obtained system satisfying the conditions 


(<?, t) = V±J(T, t ) + 0((7 °°), a -> ±CX), 


J m,s,p,q\ 




(2.13) 
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where 


V 


± 

m, s,p,q 


>P) 


t s PfL m (a), s f: m, p = q = 0, 
0 , otherwise, 


(2.14) 


the symbol 0(c r~°°) denotes a function of order 0(cr~ N ) for any N > 0. 

Denote by £ the class of C'’“-smooth functions v(a,t ) for a G M, Ti ^ t ^ T 2 . By 97t + 
we denote the set of functions from £ such that there hold the inequalities 


d t+j v(a, t) 
da i dti 


< exp(—ycr), 


Vi, j. 


By 971 we denote an analogous set of function, for which there hold the inequalities 


d l+j v(a, t) 
da i dP 


< M i>:j exp (ycr), 


Vi, j- 


Further, we need the following statement. 

Lemma 1. Let 

P- G £, P + G £, F G C, 
F-L 3 P~e9n-, F-L 3 P + eWl + . 
Then for the existence of a solution of the problem 

l 3 v = f, v-p ± em ± 


it is necessary and sufficient the fulfillment of the condition 

{ F( p+ - p~) + k- v(i-o)](p + - />-)} = 

0 oo 

= J [F(a, t) — L 3 P~(a, t)\da + J[F(a,t) - L 3 P + (a,t)\da. 

—oo 0 


(2.15) 


The proof of this lemma can be found in [3J chapter VI], 

Theorem 2. For cel and t > 0 there exist solutions of equations (12.81) - (12.1011 such 
that v 0 ~ i'o G fOF and - V^ n -m,p,q e w/iere are functions 

(12.140 . Under the condition that all v m q n r_ m > p i q r for n! < n are already determined, each 
function v m:n _ m: p :q is determined uniquely up to a term x m ,n-m,p,q v o( <7 )-> where Xm,n-m,p, q is 
an arbitrary constant. 

Proof. Formula (12.71) gives the function uo(cr) satisfying equation (12.8)) and condition 
Vq — Vq G DJLV Thus, problem (]2.9j) . (12. 13jl has the solution y 10iPi9 (cr, t) = >Cip,p,q{t)v'o(a). 
Proceeding by induction, suppose that solutions of problems (12.10)) . ()2. 13j) for n < k are 
constructed so that for n = k the problem has a solution 


V m,k—m,p,q{ (T i ^) "b ^m,k-m,p,q{t)Vo(a) . 


Let us find a function >Cm,k-m,p,q(t) for which problem 02 .10)) , ()2.13j) has a solution for n = 
k + 1. Apply Lemma 1, whose conditions are fulfilled for 


F(a,t) = 


dVm,k—m,p,q dQm,k+l—m,p,q ,, ! 


dt 


+ 


, P ± (a,t) = V± k+1 _ mpq (a,t), 


da 






























to the function v mi k+i -m,p,g(<T t). According to Lemma 1, for the solvability of problems (12.101) . 
(12.131) for n — k + 1 it is necessary and sufficient the fulfillment of conditions (j2. 15j) . Substi¬ 
tuting functions v mjn _ m:P:q (a,t) for n < k and 

v m,k—m,p,qi cr j t) V m,k—m,p,q( < 7 i^) T ^ i m,k-m,p,q{t) v O^) 

into this equality, we see that the function x mj k-m,p,q{t) enters in the form 


o 



— OO 


d(^(m,k—m,p,q(t')'U o(v)) 
dt 


dcr + 


f 9^3i rnj k—m,p,q{t')v q ( o "))^ 

J dt 


o 


dx, 


UU 


rn,k—m,p,q 


(t) 


dt 


Thus, condition (I2.15|) has the form 


/ 

^ m,k—m,p,q 


(t) = fit ), 


where f(t) is a known function. This implies that a solution v mj k-m,p,qio' 1 t) is determined 
up to a term x m ,k-m,p,q v b( (T ) an d problem (12. lGj) . (j2. 13j) is solvable for n — k +1. Theorem 2 
is proved. 


3 Dispersive compression wave 


Dispersive compression waves are studied in plasma, fluids, and optics [14, [15], lbiUZJ- hr 
particular, studying properties of solutions of the Korteweg-de Vries equation 


du du d 3 u 

dt +U fa +£ df? = ’ 


t ^ 0, e > 0, 


(3.1) 


is of indisputable interest for describing nonlinear wave phenomena. For steplike initial 
data, an asymptotic solution of the KdV equation (13. lj) can be found by the inverse spectral 
transform method [IS]. Under certain restrictions on the initial function, the asymptotic 
behavior can be studied by the Whitham method, as in Haem- The long-time asymptotic 
solution with step-like initial data is also analyzed in EBEg. 

Open questions concerning the behavior of solutions are still the subject of attention for 
modern researches [23]. Some mathematical results about the solution of the problem in 
various cases can be found in [24] and [25] . 

Here, we consider results of paper [25] for the Cauchy problem with the initial condition 


u{x, 0, e, p) = A(xp ), t — 0, i6 K, p > 0, 


(3,2) 


introduced in the previous section. We assume the fulfil lm ent of the condition 

p = ——y 0 . 

V £ 


In paper [26], an asymptotic solution to problem (13.ID (I3.2D is constructed using renor¬ 
malization EH EH E9]. Let us pass to the inner variables 


x 


= \feij, t = y/e 6, 
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since this allows one to take into account all terms in equation (13.ID . As a first approximation, 
we take the solution of the equation 


dZ dZ d 3 Z 

dO dr] dr] 3 


(3.3) 


with the initial condition 


z (v, o) 


A 0 , v < 0, 
A+, 7] > 0, 


(3.4) 


where Aq = lim A(s). As a model of collisionless shock waves, problem (j3.3jl fl3.4fl was 

s —>-dzoo 

studied by A.V. Gurevich and L.P. Pitaevskii in [3Uj . 

Let us construct the expansion of the solution in the following form: 


u(x,t,£,p ) = Z(j],6) + pW(r],0,p) + O(//*), a > 0, 


(3.5) 


where the addend pW (■ rj , 6 , p) must eliminate the singularity of Z at the initial moment of 
time. Then the function W satisfies the linear equation 


dW d(ZW) d 3 W n 

-1— -- H-= 0. 

dO dr] dr] 3 

Differentiating equation (13. 3 p with respect to 7], we find that the expression 

1 dZ{r],0) 


(3.6) 


G(v,0) = 


- A 0 - 


dr] 


satishes equation (13.6D . Moreover, G is the Green function, because 

OO CO 

1 


lim / G(j],0)f(j]) dr] — 


e^+o 


A 0 + -A 0 - 


z(v,Q)f'(v) d o = m 


for any smooth function / with compact support, thus G(r], 0 ) = 5(7]). 

Let us choose the solution W in the form the convolution with the Green function G 
so that the asymptotic approximation would satisfy the initial condition (13.21) . As a result, 
expansion (j3.5j) becomes 


u(x, t, e, p) = U Q (x, t, e, p) + 0(p a ), 


where 

CO 

U 0 (x, t, £, p) = Z(p, 0) + A+ 3 A _ J [A(s) - Z(s, 0 )] ds. 

—OO 

Integrating by parts, we obtain the asymptotic approximation 

OO 

u(x, t, £, p) S 3 Uq(x, t, £, p) = A + ^ A - J z A'(s) ds. 

—OO 

Constructing complete asymptotic expansions of the solution near the singular point by 
the standard matching method may be connected with serious difficulties. In fact, it is 
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necessary to solve the scattering problem for a recurrence system of partial differential equa¬ 
tions with variable coefficients [3]. In addition, the investigation of the shock wave generated 
by gradient catastrophe shows that the asymptotics of the solution in a neighborhood of a 
singular point may have a multiscale structure j4] as we have seen in Section 1. 

The renormalization approach allows one to construct a uniformly suitable asymptotics 
in the whole domain of independent variables avoiding difficulties arising from the matching 
procedure. 

In particular, for A$ = 0 and = a > 0 the following formula was obtained in 
paper [25]: 


u(x, t, £, p ) ~ 2A 


x + at 
P 


A 


x — 2at/3\ 


P 


) 


2/3 

— J A' ~ at V ^j 

-l 


2 dn 2 


/ a 3 / 2 tcu(y) 

V ^ 



+ ° 2 (y) 


dy, 


where dn(u, a) is the elliptic Jacobi function 


dn (u, m ) = \J\ 


msm ip, 


<p(u) 


dv 


u = 




m sm" v 


u(y) = 


1 

V6 


| y- \ [i + ° 2 {y)\ |, 


2<t 2 (1 — cr 2 )/l (cr) 
E(a) - (1 - a 2 )K(a) 


K(a) and E(a) are complete elliptic integrals of hrst and second kind: 


3 y, 


K(a) 


tt/2 



dv 

2 * 2 

— a z sm v 


E(&) — 



— a 2 sin 2 v dv. 
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